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We analyze the coherence properties of a cold or a thermal neutron by utilizing the 
Wigner quasidistribution function. We look in particular at a recent experiment 
performed by Badurek et ai, in which a polarized neutron crosses a magnetic field 
that is orthogonal to its spin, producing highly non-classical states. The quantal 
coherence is extremely sensitive to the field fiuctuation at high neutron momenta. 
A "decoherence parameter" is introduced in order to get quantitative estimates of 
the losses of coherence. 

PACS: 03.65.Bz; 03.75.Be; 03.75.Dg 

1 Introduction 

Highly non-classical, Schrodinger-cat-like neutron states can be produced by coherently 
superposing different spin states in an interferometer and with neutron spin echo [1, 2]. 
We analyze here an interesting recent experiment [3] in which a polarized neutron 
crosses a magnetic field that is orthogonal to its spin, producing Schrodinger-cat-like 
states. Our main purpose is to investigate the decoherence effects that arise when the 
fluctuations of the magnetic field are considered. 

2 Squeezing and squashing 

Let us start by looking at the coherence properties of a neutron wave packet and con- 
centrate our attention on the losses of coherence provoked by a fluctuating magnetic 
field. To this end, we introduce the Wigner quasidistribution function 

W{x, k)^l-Jd^ e-^'^H {x + r - I) ' (1) 
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where x is position, p — Tik momentum and tl) the wave function of the neutron in the 
apparatus. The Wigner function is normahzed to one and its marginals represent the 
position and momentum probabihty distributions 



j dx dk W{x, k) = 1; P{x) = j dk W{x, k), P{k) = j dx W{x, k). 



(2) 



We shall work in one dimension. We assume that the neutron wave function is well 
approximated by a Gaussian 
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where S is the spatial spread of the wave packet, SkS = 1/2, xq is the initial average 
position of the neutron and po = hko its average momentum. The two functions above 
are related by a Fourier transformation and are both normalized to one. Normalization 
will play an important role in our analysis and will never be neglected. 
The Wigner fmiction for the state (3)- (4) is readily calculated 
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and turns out to be a positive function. In the language of quantum optics [4], we shall 

say that the neutron is prepared in a coherent state if (5 = (5fe = 1/ \/2 and in a squeezed 
state if (5 7^ (5fe. An illustrative example is given in Figure 1. 

Consider now a polarized neutron that crosses a constant magnetic field, parallel to 
its spin, of intensity B and contained in a region of length L. Since the total energy is 
conserved, the kinetic energy of the neutron in the field changes by AE — fxB > 0, where 
— /z is the neutron magnetic moment. This implies a change in average momentum Ak = 
m,fiB /h^kg and an additional shift of the neutron phase proportional to A = LAk/k^. 
The resulting effect on the Wigner function is W{x, k) — > W{x — A,k). 

Assume now that the intensity of the i3-field fluctuates around its average Bq ac- 
cording to a Gaussian law. This fluctuation is reflected in a fluctuation of the quantity 
A according to the distribution law 
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where a is the standard deviation. The ratio a/ Aq is simply equal to the ratio SB/Bq, 
5B being the standard deviation of the fluctuating B-field. The average Wigner func- 
tion, when the neutron has crossed the whole B region of lenght L, represents a 
"squashed" state, that has partially lost its quantum coherence: 
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Fig. 1. Above, Wigner function (5), for xo = 0, feo = 1.7- lO^^m^^; x is in units 10~^°m 
and p — k in units lO^^m"^. Prom left to right: 5 = l/\/2 (coherent state), 5 = 1 (squeezed 
state) and S = \f2 (a more squeezed state); the uncertainty principle always reads 5^5 = 1/2 
(minimum uncertainty states). Below, Wigner function in (7) for the same values of xo,feo 
and Ao = 0. From left to right: cr = 0, l/-s/2, •\/3/2; the uncertainty principle yields (5^5' = 
1/2, 1/a/2, 1 (squashing). 



This function is represented in Figure 1 for Ag = (vanishing average magnetic field) 
and increasing values of u. The above Wigner function can be calculated explicitly, but 
its expression is a bit cumbersome; however, its marginals (2) are simple: 
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Notice that the momentum distribution (9) is unaltered [|^(/c)P in (4)]: obviously, 
the energy of the neutron does not change. Observe the additional spread in position 

8' = {S^ + (t2)1/2 and notice that the Wigner function and its marginals are always 
normalized to one. The uncertainty principle yields SkS' = ^y/l + a^j5^ > 1/2. 



3 Schrodinger-cat states in a fluctuating magnetic field 

Let us now look in more detail at the experiment [3] . A polarized {+y) neutron enters a 
magnetic field, perpendicular to its spin, of intensity Bq = 0.28mT, confined in a region 

of length L = 57cm. Due to Zccman splitting, the two neutron spin states travel with 
different speeds in the field. The average neutron wavenumber is ko = 1.7 • 10^*^m^^ 
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and its coherence length (defined by a chopper) is 5 = 1.1 • 10 



By traveUing 



in the magnetic field, the two neutron spin states are separated by a distance Aq = 

m, one order of magnitude larger than S (notice the factor 2, 
absent in the definition of the previous section) . Observe that the neutron wave packet 
itself has a natural spread dt = {6^ + {ht/2mSy y^^ ~ 15cm (due to its free evolution 
for a time t ~ mL/Tiko); however, we shall neglect this additional effect, because it is 
irrelevant for the loss of quantum coherence. 

After the neutron has crossed the B-field only the +y spin-component is observed 
and its Wigner function is readily computed 



W{x, k) = — exp[-252(fc - kof] 
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(10) 

Notice that for A = (no _B-field) one obtains (5). Our interest is to investigate the 
loss of quantum coherence if the intensity of the S-field fluctuates, like in the previous 
section, yielding a random shift according to the law (6). In such a case, the average 
Wigner function reads 

Wn,{x, k)= [ dA w{A) W{x, k) = exp[-25'^{k - kof] 
J 47r 
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and the momentum distribution function yields 



P{k) = d — exp[-25'{k-koy] 
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(11) 
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Notice also that, since only the -l-y-component of the neutron spin is observed, the 
normalization reads 



N = j dxdkW^{x,k) 
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Obviously, N = 1 when no magnetic field is present (a = Aq = 0). The Wigner 

fimction (11) is plotted in Figure 2 for some values of a. The off-diagonal part of the 
Wigner function ("trustee" of the interference effects) is very fragile at high values of 
momentum. This was already stressed in [2, 3] and is apparent in the structure of the 
marginal distribution (12): the term exp(— o-2fc2/2) strongly suppresses the interference 
effects at high A:'s. 
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Fig. 2. Wigricr function (11) for Ao — 0, fco = 1.7- lO^^ni^^, 5 — 1.1 ■ 10^^'^m; x is measured in 
units 10~^'^m and p = fc in units lO^'^m"^. From above left to bottom right, a = 0, 0.5, 1.0, 1.5- 
10~^°m. Notice the asymmetry of the "wiggles" around ko and their fragility at high values of 
momentum. Observe also the slight "squashing" of the Gaussian components at high values 
of a. 



4 Decoherence parameter 

One can give a quantitative estimate of the loss of quantum coherence by introducing 
a "decoherence parameter," in the same spirit of Refs. [5]. To this end, remember that 
the Wigner function can be expressed in terms of the density matrix p as 

Wix, k) = ^jd^ e-^'Hx + mp\x - ^/2), (14) 

and that Tr(/9^) = Trp = 1 for a pure state, while Tr(p^) < Trp = 1 for a mixture. 
Define therefore the decoherence parameter 

Tr(p2) 2nJdxdkW^{x,kf 

sia) = 1 — — — -7T = I 5". (151 

(Trp)2 {JdxdkW^{x,k))^ 
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Fig. 3. Decohcrencc parameter. Left: e as a function of S and a (both in units 10~^°m). 
Notice the peculiar behavior when 6 > 3 and 1 < a < 2. Right: e vs cr (in 10~^°m) for 
S = 1.1 ■ 10~'^°m (experimental value in [3]). 



This quantity is expected to vanish for cr = (no fluctuation of the B-field and quantum 

coherence perfectly preserved) and to become unity when cr ^ oo (large fluctuations of 
the 5-field and quantum coherence completely lost). Figure 3 conflrms these expecta- 
tions, that can also be proven anahticaUy from (11). In Ref. [3], ^ = 1.1 • 10~^°m and 
a is (presumably) very small, being the intensity of the B fleld controlled with high 
accuracy. It is remarkable that the decoherence parameter is not a monotonic function 
of the noise a, when 6>3 - 10~^°m and 1 • 10~^°m< cr < 2 • 10~^°m. This may be due 
to our very deflnition (15) or to some physical effect we do not yet understand. 
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